Energy— Momentum Tensor for the Electromagnetic Field in a Dielectric 



Michael E. Crenshaw and Thomas B. Bahder 
RDMR-WSS, Aviation and Missile RDEC, US Army RDECOM, Redstone Arsenal, AL 35898, USA 

(Dated: September 28, 2010) 

The total momentum of a thermodynamically closed system is unique, as is the total energy. 
Nevertheless, there is continuing confusion concerning the correct form of the momentum and the 
energy-momentum tensor for an electromagnetic field interacting with a linear dielectric medium. 
Rather than construct a total momentum from the Abraham momentum or the Minkowski mo- 
mentum, we define a thermodynamically closed system consisting of a propagating electromagnetic 
field and a negligibly refiecting dielectric and we identify the Gordon momentum as the conserved 
total momentum by the fact that it is invariant in time. In the formalism of classical continuum 
electrodynamics, the Gordon momentum is therefore the unique representation of the total momen- 
tum in terms of the macroscopic electromagnetic fields and the macroscopic refractive index that 
characterizes the material. We also construct continuity equations for the energy and the Gordon 
momentum, noting that a time variable transformation is necessary to write the continuity equa- 
tions in terms of the densities of conserved quantities. Finally, we use the continuity equations 
and the time-coordinate transformation to construct an array that has the properties of a traceless, 
symmetric energy-momentum tensor. 



I. INTRODUCTION 

The energy-momentum tensor is a concise way to 
represent the conservation properties of an unimpeded 
flow fleld. For most types of simple flows, the energy- 
momentum tensor is well-defined, with the notable ex- 
ception of the electromagnetic field in a linear dielectric 
material. The Abraham-Minkoswski controversy [H-ll] 
for the momentum of electromagnetic fields in a dielec- 
tric began with the derivation of the energy-momentum 
four-tensor by Minkowski [l^l- Noting that the unsym- 
metrical Minkowski tensor does not support conserva- 
tion of angular momentum, Abraham 11] proposed an 
energy-momentum tensor that was symmetric, but at 
the expense of a new phcnomcnological force. In order 
to address this constraint and additional issues, Einstein 
and Laub Q , Nelson Q , and others proposed variants of 
the energy-momentum tensor. 

The crux of the Abraham-Minkowski controversy is 
whether the electromagnetic momentum density in a di- 
electric is of the Minkowski form 

gM-^(DxB) (1) 

or the Abraham form 

g^ = i(ExH). (2) 

c 

Experimental efforts to resolve the theoretical impasse 
have not been definitive. While some experiments fa- 
vor the Abraham formula, other experiments support 
Minkowski's version. Brevick's [l2| analysis of exper- 
iments performed by Jones and Richards [l^, Ashkin 
and Dziedzic [3], and others showed that the allocation 
of momentum between the field and material was the 
determining factor in whether a particular experimental 
result was described by the Abraham or Minkowski form 
of electromagnetic momentum. Following Brevick [l^ . 



the formula for the field momentum has been shown re- 
peatedly to be arbitrary such that any of the formulas 
for the field momentum can be combined with an appro- 
priate momentum for the material to produce the same 
total momentum 

In 1973, Gordon fl*] constructed the total momentum 
from a microscopic model in which the electromagnetic 
field component of the total momentum is said to be the 
Abraham momentum and the dielectric is treated as a di- 
lute collection of electric dipoles with center-of-mass mo- 
tion in the direction of propagation of the field. Gordon 
1] discusses the empirical and experimental validation of 
the total momentum density 

Tl 

gG = -(E X B) (3) 

c 

and shows that the density go, integrated over a volume 
containing the entire field, is invariant in time. Gordon 
concludes that the Abraham momentum density repre- 
sents the true momentum density of the electromagnetic 
field and the Minkowski momentum density includes a 
pseudomomentum. 

In this article, we investigate conservation of energy 
and momentum in a thermodynamically closed system 
consisting of the macroscopic electromagnetic field and 
a negligibly reflecting linear dielectric (such as a dielec- 
tric with an anti- reflection coating). Here, we identify 
the Gordon momentum [ij as the total momentum by 
the fact that it is invariant in time and therefore a con- 
served quantity in an isolated system. We find that it is 
not necessary to decompose the Gordon momentum into 
the sum of a field momentum, such as the Abraham or 
Minkowski momentum, and a material momentum, such 
as the canonical or kinetic momentum [13]. Instead, 
we work with the total energy and total momentum and 
derive continuity equations in terms of the densities of 
these conserved quantities. The continuity equations are 
then used to construct an array that has the properties 
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of a traceless, symmetric energy-momentum tensor, but 
in a coordinate system with time-like coordinate ct/n. 



II. ENERGY-MOMENTUM TENSOR OF 
NONINTERACTING PARTICLES 

In the continuum hmit, the density of any property of 
identical noninteracting particles can be treated as the 
number density multiplied by the amount of the property 
that can be attributed to each particle. The continuity 
equation corresponding to a specific property, such as 
mass, charge, or energy, is then obtained by substitution 
of the specific property density for a placeholder number 
density. For an infinitesimal element of volume in an 
inviscid sourceless flow, the continuity equation 



ds_ 

dt 



V • su = 



(4) 



is derived by applying the divergence theorem to a Taylor 
series expansion of the property density field s and the 
vector velocity field u = {u^ , ,u^ ) of the flow [l^. The 
continuity equation reflects the conservation of a contin- 
uous scalar property in a flow in terms of the equality of 
the net rate of flux out of the volume and the time rate 
of change of the property density s inside the volume. 
Depending on the context, the three- vector 



g = su 



(5) 



is known as the momentum density, the flux, or the cur- 
rent density of the property. Specifically, g corresponds 
to the linear momentum density if s is a mass density 
and corresponds to the charge current density if s is the 
electric charge density. 

Some conserved properties, such as momentum, are 
vectors. For a flow, the density of a conserved vector 
property can be represented as s = [s^ , , s^). Applying 
the scalar formalism to the three orthogonal components 
of the property density vector yields scalar continuity 
equations 



dsy 

'dt 
'dT 



V-s^u = 



V • s^u = 0. 



(6a) 



(6b) 



(6c) 



At this point, we adopt a four-dimensional notation 
where repeated indices are summed. We take Roman in- 
dices to run from I to 3 and we identify the coordinates 
with the Cartesian coordinates, such that = x, 
x'^ — y, and x^ = z. Greek indices run from to 3 
and is identified with the time-like coordinate ct. The 
Minkowski space-time metric is diag{—l, 1, f , 1). Finally, 



partial differentiation with respect to the indexed coor- 
dinates is represented by da = d/{dx°'). 

The four continuity equations Q and ^ can be con- 
cisely represented by 



Q/3 _ 



0, 



where 



c s 


cu^s 


cuy s 


cu^s 


cs^ 




uys"" 


u^s^ 


csy 


u^sy 


uysy 


u'^sy 






uys'' 


u^s^ 



(7) 



(8) 



Consider a transformation to a new set of coordinates 
= fi{x^,x^,x^,x^). In order for this matrix to trans- 
form as a tensor, s^, s^, and must be expressible in 
terms of s and components of u. For a closed system, 
conservation of angular momentum requires G"^ to be 
symmetric [loj 



G" 



(9) 



Therefore, the vector property s must represent the flux 
of the conserved scalar quantity s such that s = su. Con- 
versely, the flux su of a scalar property of the particles 
is a conserved vector property of the flow. The covariant 
form of the continuity equation is the four-divergence, 
Eq. ([7]), of the energy-momentum four-tensor 



c s csu 
csu^ su^u" 



csuy 

csu^ 



suyu^ 

su'^u^ 



csu" 
su^u" 

suyuy 



csu 
su^u^ 

suyu^ 

su^u^ 



(10) 



The continuity equation ([7]) with the four-tensor ([TO]) is 
valid for any conserved extensive quantity in a simple 
flow, not just for the mass of a fluid. 

The energy-momentum tensor given in Eq. (jlOp has 
some essential properties. First, the four-divergence of 
each row vector, Eq. Q, is a continuity law correspond- 
ing to the conservation of the property represented by the 
property density in the first element of that row. Second, 
the tensor is symmetric, corresponding to the absence of 
unbalanced shear forces and conservation of angular mo- 
mentum in a closed system p^ . Third, as a consequence 
of diagonal symmetry, 



daG' 



0, 



(11) 



the four-divergence of each column vector is a continuity 
law corresponding to the conservation of the property 
represented by the property density in the first element 
of that column. It should also be noted that this simple 
energy-momentum tensor is based on the properties of 
an unimpeded flow. If the flow is redirected by impact 
with a macroscopic object then one is obligated to include 
the equations of motion of the object or the forces of 
restraint. 
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III. THE ABRAHAM AND MINKOWSKI 
ENERGY-MOMENTUM TENSORS 

The Abraham and Minkowski energy-momentum ten- 
sors are examples of a number of different tensors that 
have been proposed for the electromagnetic field in a di- 
electric [15]. The Minkowski tensor can be constructed 
from the continuity equations for energy flux and momen- 
tum flux by the same procedure that was used to con- 
struct the array ([S]). The Minkowski tensor is not sym- 
metric and therefore violates angular momentum conser- 
vation if it is the total energy-momentum tensor of a 
closed system. Consequently, the Minkowski tensor is 
considered to be a representation of the energy and mo- 
mentum of a component of the system [l5| . In this sec- 
tion, we outline the construction of the Minkowski ten- 
sor and describe how the procedure is modified to obtain 
the Abraham tensor. Neither the Minkowski tensor, nor 
the Abraham tensor, satisfy the requirements of a to- 
tal energy-momentum tensor. The Minkowski tensor is 
not symmetric while the Abraham tensor contains a phe- 
nomenological volume force. 

The macroscopic Maxwell equations of continuum elec- 
trodynamics are the basis for deriving continuity equa- 
tions for electric and magnetic fields in a dielectric. For 
a dielectric with no free charges in a regime of negligible 
absorption and dispersion, the Maxwell equations may 
be written as 



V X E = - 



dB 



V X B 



d{ct) 



V-B = 



-E = 



(12a) 



(12b) 



(12c) 



(12d) 



in Heaviside-Lorentz units. The electric and magnetic 
fields can be defined in terms of the vector potential A 

as 



E 



dA 

diet) 



B = V X A 



(13a) 



(13b) 



for transverse fields in the Coulomb gauge. 

The macroscopic Maxwell equations are the axioms of 
classical continuum electrodynamics. Poynting's theo- 



d 



diet) 



2 ^ 



B^ 



V • (E X B) = 0, (14) 



can be derived by multiplying the Faraday law (I12ap by 
B and adding it to the Maxwell- Ampere law (jl2bp mul- 
tiplied by E. Poynting's theorem can also be derived 
by substituting the Maxwell equations into the temporal 
derivative of the energy density 



(l/2)(n^E'' + B 



(15) 



using a vector triple-product identity. Poynting's theo- 
rem is a continuity equation for the Poynting energy-flux 
vector 



Sp = c(E X B) = (s^,sp,4) 



(16) 



The theorem 
d 



diet) 



{n^E X B) = -B X (V X B) + B(V • B) 



n^E X (V X E) +E(V -n^E) 



(17) 



is derived, in a manner similar to Poynting's theorem, 
by substituting Maxwell's equations into the temporal 
derivative of the Minkowski momentum density, gM = 
(n2/c)(E X B). The right-hand side of Eq. (H?]) can be 
recast, approximately, as the negative of the divergence 
of the Maxwell stress tensor ,20] with components 



-n^EiEj - B,Bj + - (n^E • E + B • B) (5^, 



(18) 

where terms involving the gradient of tt? have been ne- 
glected. Then the temporal derivative of the Minkowski 
momentum density, Eq. (jl7l) , can be expressed using the 
vector divergence operator, V •, as 



d 



The array 



d{ct) 



pa/3 _ 
'-M — 



(n2ExB)-^V-W = 0. 



Pe s\i/c Sp/c Sp/c" 

cglj M^" W^^ W^^ 

Cglt ^22 1^23 

cgli VF^l 1^32 |y33_ 



(19) 



(20) 



known as the Minkowski energy-momentum tensor, is 
constructed from continuity equations (|14p and (flQl) . Us- 
ing the summation convention, we can write Poynting's 
theorem in Eq. as 







for continuity of the energy flux and Eq. (fT9| as 



(21) 



(22) 



for continuity of the momentum flux. Then, each row of 
Tm corresponds to a four-divergence 



0. 



(23) 
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Because the array T^f in Eq. (|20l) is not symmetric, 
Abraham proposed the energy-momentum tensor 



T 



.Q/3 



Pe 

c9a 
c9a 



Sp/C 


4/c 












VF22 


Vl/23 




^^32 


Vl/33 



(24) 



where gA is given by Eq. ([2]). The four divergence of this 
tensor is 



-f 



(25) 



where is the Abraham force. The Cartesian compo- 
nents of the Abraham force 



d 



diet) 



((^2 - 1)E X B) 



(26) 



are obtained by substituting the Abraham momentum 
density into the continuity equation (jl9p and the time- 
hke coordinate component is /'^ = 0. 

It has been widely reported in the hterature that nei- 
ther Tm nor Ta is the total energy-momentum tensor 
[is'l. Instead, they are to be considered two of many ar- 
bitrary forms of the electromagnetic part, Tfld, of a total 
energy momentum tensor 



matl 



composed of energy-momentum tensors for the field and 
material subspaces, however those subspaces are defined. 



IV. TOTAL MOMENTUM 

A well-defined quantity for the total momentum is 
only derivable from the macroscopic Maxwell equations 
by imposing an additional condition: The total momen- 
tum of an isolated system must be constant in time. 
In this section, we identify the unique total momentum 
using this constraint. We consider the case of a quasi- 
monochromatic electromagnetic field, in the plane-wave 
limit, entering a linear medium from vacuum at normal 
incidence. The medium is taken to be a simple linear 
dielectric in the regime of negligible dispersion and neg- 
ligible absorption. In the absence of reflection, there is 
no momentum given to the dielectric slab and it remains 
stationary. We adopt this case of a stationary dielectric 
in which reflections can be neglected by assuming that 
an antireflection coating has been applied to the dielec- 
tric or that the refractive index of the dielectric is only 
slightly greater than unity. 

Gordon [ij used a microscopic model of the dielectric 
as a vapor of weakly polarizable atoms and derived the 
material momentum as the continuum average of the me- 
chanical momentum of the atoms. Assuming a rarefied 
vapor of atoms in order for reflections to be negligible, 
Gordon obtained the total momentum 

Gg= / dvgG^ I di'-(ExB) (28) 
Jv Jv c 



by adding the material momentum to the Abraham 
momentum. Here the integration is over all three- 
dimensional space. In continuum electrodynamics, the 
electrodynamic properties of a material are character- 
ized only by a macroscopic refractive index. Therefore, 
the microscopic origin of the Gordon momentum ()28|) is 
of no consequence in the formalism of continuum electro- 
dynamics. Following Gordon, we will demonstrate that 
the momentum (l28t is invariant in time and, because it 
is conserved, can be identified as the total momentum. 
Because the Gordon momentum depends on the material 
only through the refractive index, it is the unique total 
momentum for all cases in which the medium behaves, 
to a good approximation, as a negligibly reflecting linear 
dielectric with refractive index n. 

Propagation of a fleld in a linear medium is governed 
by the wave equation. 



n2 d'^A 
"^"5*2" 



0, 



(29) 



written in terms of the vector potential A, where B — 
V X A. For quasi-monochromatic plane waves, it is con- 
venient to write the vector potential in terms of a slowly 
varying envelope function, A(z,t), a rapidly varying car- 
rier, and a unit vector, et, in the direction of propagation 



(27) as A = A{z, t)e '^^•'et. For a plane-wave entering a 



dielectric at normal incidence, reflections are negligible if 
6n — n — I is small. In this limit, there is no momentum 
given to the bulk material, which remains stationary, and 
we can apply the Fresnel relation 



A 

A,. 



= [! + - 
n+1 \ 2 



(30) 



where Ai is the incident amplitude and At is the trans- 
mitted amplitude. Comparing Eq. (j30l) with a series ex- 
pansion of l/Vn in the limit of small 6n, we flnd that 
the vector potential amplitude inside the dielectric. At, 
is reduced by a factor of ^/n from the incident amplitude 
A,. 

For continuous plane waves, and approximately for 
slowly varying waves, the relation between the ampli- 
tudes of the flelds simplifles to 



|B|=n|E| = — lAI 



(31) 



Then, the electromagnetic energy density, pe 
(l/2)(n2E2 -I- B2), can be written as 



Pe = 



(32) 



Applying the Fresnel amplitudes to relate the fields in- 
side and outside the medium. At = Ai/ y/n, we find that 
the energy density inside the material is a factor of n 
larger than the energy density of the same field in the 
vacuum. The region occupied by the field in the mate- 
rial is compressed spatially by a factor of n due to the 
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reduced speed of light within the medium, such that the 
total energy 



U : 



dvpe 



dv- 



(33) 



is conserved. Numerical solutions of the wave equation 
for a field entering a linear material through a gradient- 
index anti-reflection coating indicate that the field in the 
material is a factor of ^/n smaller and a factor of n nar- 
rower than the field in the vacuum, independent of the 
magnitude of n, as long as reflections are suppressed [2l| . 

Having demonstrated the conservation properties of 
the electromagnetic energy, we demonstrate the conser- 
vation properties of the Gordon electromagnetic momen- 
tum by a similar procedure. The Gordon momentum 
[H is obtained by integrating the momentum density ([3]) 
over all three-dimensional space. Comparing the Gordon 
momentum, expressed in terms of the envelope functions 



Gr 



dv- (E X B) 

V c 



dv- 



efc, 



(34) 



with the total energy in Eq. we see that conserva- 
tion of total energy implies conservation of the Gordon 
momentum. The total momentum of a closed system is 
unique and the Gordon form of total momentum is con- 
served. Therefore, the Gordon momentum can be iden- 
tified as the total momentum of the thermodynamically 
closed system. Because neither the Minkowski momen- 
tum nor the Abraham momentum is the total momen- 
tum, neither is conserved in a thermodynamically closed 
system. 

We point out that, in the macroscopic limit in which 
the dielectric is described by a refractive index n, the 
question of what portion of the energy given by Eq. ([33|) 
resides in the field or dielectric is improperly posed. Com- 
paring Eqs. (p3)) and (p4)) . we see that the same holds true 
for the apportionment of the momentum into field and 
dielectric components. 



V. TOTAL ENERGY-MOMENTUM TENSOR 

In the previous section, we identified the unique to- 
tal momentum for the system of an electromagnetic field 
in a dielectric for the case of negligible reflections. In 
this section, we construct the corresponding total energy- 
momentum tensor from continuity equations for the en- 
ergy and momentum. 

A continuity equation is a differential form of a con- 
servation law applied to an element of volume in a con- 
tinuous flow. In Sec. II we showed that the energy- 
momentum tensor for dust is constructed from conti- 
nuity equations in which the differential operators act 
on the densities of conserved quantities. The operand 
of the time derivative in the continuity Eq. (fT9|) is the 
Minkowski momentum density. Because the volume in- 
tegral of the Minkowski momentum density is not a con- 
served vector quantity, we do not consider Eq. (|19p to be 



a suitable continuity equation with which to construct an 
energy-momentum tensor. Instead, we write Eq. (1191) as 



n d 
c9t 



cgG +V- W = 0, 



(35) 



in terms of the Gordon momentum density in Eq. ([3]) 
whose volume integral is a conserved vector quantity. 
Equation p5|) provides three continuity equations for 
our energy-momentum tensor. The additional continu- 
ity equation is obtained by writing Poynting's theorem 
in ^ as 

-^ + V-[7i(ExB)] = — •n(ExB) (36) 
c ot n 

using the densities of conserved quantities. Again, we are 
considering the case of a closed system in which there are 
no reflections. As this is assumed to be accomplished by 
a gradient-index antireflection coating, we can drop the 
term containing Vn and write Poynting's theorem in Eq. 
(I36t as 



c dt 



V • [n(E X B)] = 0. 



(37) 



In Maxwell's equations, expressed in terms of 3- 
vectors, time is not a coordinate. We make a simple 
change of time variable to t = t/n and write the conti- 
nuity theorems (|55|) and ([57)) as 



1 d 

C OT 



'•W==0 



1 dpe 

c dr 



V • CgG = 0, 



(38a) 



(38b) 



where the Gordon momentum density gc is given in Eq. 
^ . However when writing Maxwell's equations as tensor 
equations, time is one of the four space-time coordinates 
and we define the time-like coordinate 



x° = cT^-. (39) 
n 

Then the four scalar continuity equations, Eqs. pSap and 
(j38b[) . can be written concisely as a single equation, as 
in Section 2, defining the operator 



da 



UxO' 



dx:dy,d2 



(40) 



and an array 



T 



Pe 

-C5g 



c5g cgl 



1^21 



c9g 

1^23 



such that 



(41) 



(42) 
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The array that appears in Eq. (j4ip has a number of 
notable properties. The array is symmetric 

and has a vanishing trace 

= 0. (44) 

The operator defined in Eq. (|40|) applied to the rows 
of the array in Eq. (I4ip generates continuity laws for 
demonstrably conserved electromagnetic energy and mo- 
mentum properties. Similarly, the operator in Eq. (j40p 
applied to the columns, 

a„T"'5 = 0, (45) 

generates the same continuity equations as a consequence 
of the symmetry of the array. These are the proper- 
ties that we associate with an energy-momentum tensor. 
Ravndal [l^j arrives at the same energy-momentum ten- 
sor as in Eq. (j4ip using symmetry arguments, but inter- 
prets the continuity equations, Eqs. (j42p . in the context 
of the Minkowski momentum, which we have shown is 
not conserved. 



VI. CONCLUDING REMARKS 



momentum tensor that is comprised of separate field and 
matter tensors — an approach that involves assumptions 
about the behavior of matter in the presence of an elec- 
tromagnetic field. In continuum electrodynamics the in- 
teraction of the field and matter is described in terms of a 
single macroscopic parameter, the refractive index n. We 
showed that the Gordon momentum is the total momen- 
tum of a thermodynamically closed system consisting of 
a quasimonochromatic field and negligibly reflecting lin- 
ear dielectric in the continuum. We derived continuity 
equations from the Maxwell equations and used a time 
variable transformation to write the continuity equations 
in terms of densities of conserved energy and Gordon mo- 
mentum quantities. When written in four-dimensional 
tensor form with time- like coordinate ct/n, the continu- 
ity equations are obtained from the four-divergence of a 
traceless, diagonally symmetric energy-momentum ten- 
sor. 

In summary, for the case of an electromagnetic field 
and negligibly reflecting dielectric, we constructed the 
total energy-momentum tensor in Eq. ()4ip from conti- 
nuity equations that were derived from the macroscopic 
Maxwell's equations. It is interesting to note that a time 
coordinate transformation was required in order to write 
the continuity equations as the four- divergence of the 
symmetric energy-momentum tensor. 



For many years, the Abraham-Minkowski contro- 
versy has been resolved by postulating a total energy- 
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